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Many practical applications exploit an external loca l magnetic field - mag netic obstacle 



- as an essential part of their construction. Recently, IVotvakov et all (|2007f ) have demon- 
strated that the flow of an electrically conducting fluid influenced by an external field 
can show several kinds of recirculation. The present paper reports a 3D numerical study 
whose some results are compared with an experiment about such a flow in a rectangular 
duct. First, we derive equations to compute analytically an external magnetic field and 
verify these equations by comparing with experimentally measured field intensity. Then, 
we study flow characteristics for different magnetic field configurations. The flow inside 
the magnetic gap is dependent mainly on the interaction parameter N, which represents 
the ratio of the Lorentz force to the inertial force. Depending on the constrainment factor 
k = My/ L y , where M y and L y are half-widths of the external magnet and duct, the flow 
can show different stationary recirculation patterns: two magnetic vortices at small k, 
a six- vortex ensemble at moderate K, and no vortices at large n. Recirculation appears 
when N is higher than a critical value N Cj7n . The driving force for the recirculation is 
the reverse electromotive force that arises to balance the reverse electrostatic field. The 
reversion of the electrostatic field is caused by a concurrence of internal and external 
vorticity correspondingly related to the internal and external slopes in the M-shaped 
velocity profile. The critical value of N c ^ m quickly grows as n increases. For the case of 
well developed recirculation, the numerical reverse velocity agrees well with that obtained 
in physical experiments. Two different magnetic systems induce the same electric field 
and stagnancy region provided these systems have the same power of recirculation given 
by the N/N c ^ m ratio. The 3D helical peculiarities of the vortices are elaborated, and an 
analogy is shown to exist between a helical motion inside the studied recirculation and 
a secondary motion in the process of the Ekman pumping. Finally, it is shown that a 
2D model fails to properly produce stable two and six-vortex structures as found in the 
3D system. Interestingly, these recirculation patterns appear only as time dependent and 
unstable transitional states before the Karman vortex street forms, when one suddenly 
applies a retarding local magnetic field on a constant flow. 



1. Introduction 

An electrically conducting fluid flow influenced by a local external magnetic field is 
of considerable fundamental and practical interest. Applied to the flow, a transverse 
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homogeneous magnetic field creates a so-called magnetic obstacle, i.e. a region where the 
flow motion is retarded by the Lorentz force. 

On the fundamental side, such a system possesses a rich variety of dynamical states. 
This can already be imagined from the fact that its behavior is characterized by two pa- 
rameters, the Reynolds number Re — u$H/v and interaction parameter N = aH Bq/ puq, 
where H, uq, B are characteristic scale, velocity and magnetic field induction, and p, u, 
a are d ensity, kinemat ic viscosity, and e lectric conductiv ity of the fluid, see e.g. ISherclifi 
( 19621) ; iRobertd (jl967t) ; iMoreaul (|l990h ; iDavidsonl (|200l[ ). Re represents the ratio of the 
inertial to viscous forces in the flow, and ./V represents the ratio of the Lorentz forces to 
the inertial forces. In ordinary hydrodynamics, such as a flow around a solid obstacle, an 
increase of the inertial force, i.e. Re, renders a system to the nonlinear dynamics charac- 
terized by a vortex motion past the obstacle, see Fig. [TJi. The additional flow parameter 
TV brings new nonl i near d egrees of freedom into the problem as was elaborated recently 
bv lVotvakov et all (|2007[) . see Fig. [TJ>. 

On the practical side, spatially localized magne tic fields p l ay an essential role in a va- 
riety of industrial applications in metallurgy, e.g. iDavidson (1999), includin g stirring of 
melts by a moving magnetic obstacle (called electromagnetic stirring, e.g. iKunstreichl 
(2003)), removing undesired turbulent fluctuations during steel cas t ing u sing steady 
magnetic obstacles (called electromagnetic brake, e.g. iTakeuchi et al. (2003)) and non- 
con tact flow me a surem ent using a magnetic obstacle (called Lorentz force velocimctry, 
e.g. Thess et al. ( 20061 )). For instance, it is important to understand whether the use- 
ful turbulence-damping effect of a magnetic brake is obliterated by excessive vorticity 
generation in the wake of the magnetic obstacle. 

As is well known in the flow past a solid obstacle, there is a stagnancy region where 
one can observe a recirculation in the appropriate Reynolds number regime, so called 
attached vortices shown in Fig. [TJi. If instead of the solid obstacle one places a magnetic 
obstacle, by means of a local external magnetic field, then there will appear electrical 
eddy currents j which induce the Lorentz force Fl = j x B. The largest retarding effect 
occurs where the transverse magnetic field B is a maximum. Therefore, past the mag- 
netic obstacle there is also a stagnancy region where a kind of a reverse flow might be 
obtained. This analogy between a solid and magnetic obstacle has been realized from the 
beginning of magne tohydrodynamics (MHD). In the earliest 2D numerical simulation, 
Gelfgat et al. ( 1978j ) have observed a kind of recirculation and marked an analogy with a 



solid body: 'the qualitative pattern of the streamlines in such a flow is similar to the situ- 
ation which arises in th e case of the flow around obj ects '. However, the specially designed 
physical experiments of lGelfgat fc Olshanskiil (|l978l ) following upon this simulation failed 
to confirm their previous numerical results: 'special attempts which we made to detect 
zones with return flow were not successful. The negative flows which occur in certain 
numerical calculations are obviously due to inaccuracies in the calculation'. The authors 
of this quotation have been correct about the inaccuracies of the refereed 2D approach 
in the sense that this approach is not suitable to describe their experiments, nevertheless 
there have been still chances to reveal a kind of recirculation in the experiments. We 
shall discuss this item later in the Section [321 



For Western readers, the term 'magnetic obstacle' has been revived in the work Cuevas et 



(|2006al ). (I n the seventies of last century, one of the authors, Yu.K., used 'magnet i c obsta - 
cle' as a working term in Riga, MHD center of the former USSR). ICuevas et al. ( 2006a ) 
have performed a 2D numerical study and described a Karman vortex street past the 
obstacle similar to those observed past a circular solid cylinder. We discuss a link between 
our 3D stationary and their 2D nonstationary results in Section 13.71 

Most recent results about the wake of a magnetic obstacle have been obtained by 
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Figure 1. Structure of the wake of a solid (a) and magnetic obstacle (b). By forming the wake, 
the solid obstacle develops attached vortices, while the magnetic obstacle develops inner (first 
pair), connecting (second) and attached vortices (third pair). 



Votvakov et al\ (|2007h . Their main result is presented in Fig. [T] where the qualitative 
structure of the wake of a magnetic obstacle is given in compari son with the wake of a 
physical obstacle. By means of 3D simulation and experiments, IVotvakov et al. ( 2007 ) 
have found that the liquid metal flow shows three different regimes: (1) no vortices, when 
viscous force prevails at small Lorentz force, (2) one pair of inner magnetic vortices 
between the magnetic poles, when Lorentz force is high and inertia small, and (3) three 
pairs, namely, magnetic as in (2), connecting and attached vortices, when Lorentz and 
inertial forces are high. The latter six-vortex ensemble is shown in Fig. 

An important factor for the flow influenced by an external magnetic field is the span- 
wise heterogeneity of the field. One can distinguish two extremal cases: (i) pointwise 
braking Lorentz force, and (ii) spanwise homogeneous braking Lorentz force. The latter 
case can be easily created, e.g. by external magnets long enough to overlap the duct, 
while the first case represents an idealization since it is impossible to have a point- 
wise external magnetic fi eld. T he first case is well studied in a 2D stratified flow, e.g. 
Voropavev &: Afanasvevl (|1994l ): it is shown that a dipolar vorticity is generated in the 
vicinity of the origin of the point braking force. Applied to MHD flows similar results 



were o btained with a 2D numerical s imulat ion of a creeping flow both in IGelfgat et 



( 1978h and recently bv lCuevas et al. ( 20066 ). The vortex dipole observed in those works 
is of the same nature as the magnetic vortices that we will discuss in the present paper. 
The second case, spanwise homogeneous magnetic fields, is of traditional MHD interests 
and well understoo d. In part i cular M-shaped profile is developed under stream wise mag- 
netic field gradien t, Sherclmj ( 1962T I. It was studied ext ensive l y in experiments dKit et al. 
( 197C), Tananaev ( 1979h ) and by numerical s imulation (Sterll ( 1990h Votvakov fc Zienicke 
( 2007 )). The mos t recent numer i cal pa per ( Votvakov fc Zienickd ( 2007 )) by comparing 
with experiment (jAndreev et al\ ( 2007 )) has established that when turbulent pulsations 
are suppressed by an external magnetic field, it is the interaction parameter N which 
governs the flow. It has been shown numerically that a spanwise homogeneous magnetic 
field is not able to reverse electric field inside the magnetic gap. As we shall prove below 
this is a necessary condition to induce recirculation between magnetic poles. 

The goals of the pr esent paper are the following: (1) report details not published in 
Votvakov et all (|2007l ) due to lack of space, (2) investigate thoroughly how a constrain- 
ment of MHD flow influences stationary vortex patterns, (3) explain the driving force for 
the recirculation, (4) find out what are the 3D peculiarities of the flow, and (5) clarify 
whether a 2D flow contains the observed stationary vortex patterns. It will be demon- 
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strated that the decisive parameter for the constrained MHD flow is a power of the 
recirculation between magnetic poles given by the N/N c ^ m ratio, where N c ^ m is a critical 
value of the interaction parameter to induce magnetic vortices at the given magnetic field 
configuration. Moreover, several successful comparisons with available experimental data 
will be given for the intensity of the magnetic field, and for a maximal stationary reverse 
flow inside the magnetic obstacle. It will be explained that the magnetic vortices firstly 
appear due to the reverse electromotive force which is induced in order to balance the 
reverse electrostatic field inside the magnetic obstacle. Also, we will discuss a 3D versus 
2D numerical approach and show that the found vortices have a 3D helical structure, 
while in a 2D model, at the given range of parameters, these vortices are not fixed by 
magnetic field and generate vortex shedding. 

The subject of the present paper has a close connection to questions of stability of MHD 
flows, however, we have not included a full bifurcation analysis of new stationary flow 
patterns. The paper sheds light on the physical factors that determine the occurrence 
of stationary recirculation, i.e. the spanwise inhomogeneity of the magnetic field and 
the necessity of a three-dimensional geometry. We consider the bifurcation and stability 
analysis as a further step, which — from the practical point of view — also involves 
additional programming work on our code to allow for the computation of Jacobi matrices 
and their eigenvalues in a high dimensional dynamical system. Therefore, besides others, 
especially the following question will remain open: are the topological changes of the flow 
patterns, which we have observed in changing the system parameters, caused by changes 
of stability, i.e. bifurcations, or not (i.e. staying on the same solution branch but the 
solution only changes topology)? This and other open questions surely deserve further 
investigation for MHD channel flow in inhomogeneous magnetic field. 

The structure of the present paper is as follows. Section [2] presents a sketch of the 
model, equations and a 3D numerical solver. As an essential part, it describes in Section 
I2.2l an analytical method to deal with magnetic field of arbitrary configuration. Section [3] 
presents the results of our numerical simulations: stationary flow patterns in the middle 
plane in Section 13.11 stability diagram in Section 13.31 mechanism for recirculation in 
Section 13.41 the 3D peculiarities of vortices in Section 13.61 as well as a relationship 
between 3D and 2D numerical methods in Section 13.71 Finally, the paper ends with a 
conclusion of the observations. 



2. Problem Definition 

2.1. Model, equations, numerical method 

A schematic of the model is presented in Fig.O a conducting fluid flows in the rectangular 
duct of dimensions Length x Width x Height = 2L X x 2L y x 2H (half-length of the duct is 
shown); x-axis corresponds to the main direction of the flow. Top, bottom and side walls 
of the duct are no slip and electrically insulating. The magnets of horizontal dimensions 
Width x Length = 2M y x 2M X are assembled symmetrically on the top and bottom walls, 
where 2h is the distance between north and south poles. The center of the magnetic gap 
is the center of the coordinate system. The constrainment factor k — M y /L y defines the 
spanwise distribution of the magnetic field; K is the varied geometric parameter in the 
present simulations. Below, we will refer to the case of k = 0.02 as a magnetic blade, 
K = 0.4 as a middle magnet, and k = 1.0 as a broad magnet. 

If nothing else is specified, throughout the paper the following geometric parameters 
have been taken: L x = 25, L y = 5, H = 1, h= 1.5, M x = 1.5, 0.02 ^ n ^ 1. Reynolds 
number, Re, and interaction parameter, N, are defined with half-height of the duct H, 
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Figure 2. Sketch of the model. Through the paper it is used the constrainment factor 
k = M y I Ly, where M y and L y are half-width of the magnet and duct, correspondingly. 



the mean flow rate Uo, and magnetic field intensity B taken at the center of the magnetic 
gap, x = y = z = 0. 

The governing equations for electrically conducting and incompressible fluid are derived 
from the Navier-Stokes equation coupled with the Maxwell equations for moving medium 
and Ohm's law. We apply the quasi-static (induction-less) approximation where it is 
assumed that an induce d magnetic field is infinitely small in comparison to the external 
magnetic field, see, e.g. iRobertsi (1967), and is therefore neglected when one calculates 
the Lorentz force, but it is not neglected when finding the electric current density j. The 
resulting equations in dimensionless form are: 



du 



+ (u- V)u = -Wp- 
j = -V0- 



1 

Re 

u x B 



Au + 7V(jxB), Vu 

v-j = o, 



0, 



(2.1) 
(2.2) 



where u denotes velocity field, B is an external magnetic field, j is electric current den- 
sity, p is pressure, 4> is electric potential. The interaction parameter N and Reynolds 
number Re, N — Ha 2 /Re, are linked by means of the Hartmann number Ha: Ha = 
HBq(o j pv) x l 2 . The Hartmann number determines the thickness of Hartmann boundary 
layers, S/H ~ Ha -1 for the flow under constant magnetic field. For a given conducting 
fluid and geometry of the duct, one varies either the flow rate velocity uq, i.e. Re, or the 
magnetic field intensity Bq, i.e. Ha. In both cases, N changes. 

At given external field B, the unknowns of the partial differential equations (|2.1I - 
12. 2p are the velocity vector field u(x,y,z), and two scalar fields: pressure p(x,y,z) and 
electric potential 4>(x, y, z). For no-slip and insulating walls, the boundary conditions are 
u\t = 0, <90/(9n|r = 0, where n is normal vector to a surface T. The outlet of the duct 
was treated as a force free (straight-out) border for the velocity. The electric potential 
at inlet and outlet boundaries was taken to be equal to zero because the inlet and outlet 
are sufficiently far from the region of magnetic field. The stationary laminar profile of an 
infinite rectangular duct known analytically in the form of a series expansion was used 
as the inlet velocity profile. Because we are interested in a stationary solution, the initial 
conditions play no role (except for the speed of convergence) . 



The 3D numerical solver has been described in details in lVotvakov fc Zienickd (|2007j) . 
It was developed from a free hydrodynamic solver created originally in the research 
group of Prof. Dr. M. Griebel (|Griebel et all (|l995| )). The solver employs the Chorin- 
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type projection algorithm and finite differences on an inhomogeneous staggered regular 
grid. Time integration is done by the explicit Adams-Bashforth method that has second 
order accuracy. Convective and diffusive terms are implemented by means of the VONOS 
(variable-order non-oscillatory scheme) scheme. The 3D Poisson equations for pressure 
and electric potential, arising at each time step, are solved by using the bi-conjugate 
gradient stabilized method (BiCGStab). 

The computational domain, |x| ^ L x , \y\ ^ L y , \z\ H, has been discre t ized b y an 
inhomogeneous regular 3D grid described in detail in Votvakov fc Zienicke (2007). To 
verify that the inlet and outlet boundaries have no influence on the presented results, 
we have carried out several simulations with double the number of grid points in the 
IE-direction and found no differences. Moreover, we have also varied the inhomogeneous 
grid resolution both in y and z-direction to be sure that Hartmann and sidewall layers 
are properly resolved. 

2.2. Fast analytical method for a proper magnetic field 

As is imposed by the electrodynamics, an external magnetic field must be divergence 
and curl-free. Although authors of previous works realized this fact, to define their fields 
they used simple mathem atical functi o ns which did not s atisfy divergence and/or curl- 
free requirements, see, e.g. ISterll(|l990l) .lAlb oussier e1 (|2004b . This fact might be explained 
by many reasons, e.g. either by the complexity of a real field or insignificance of effects 
appearing due to an inaccuracy of the field definition. Therefore we feel that there are 
insufficient correct and yet simple methods to define an external magnetic field of arbi- 
trary configurations. To fill this gap we explain below a simple physical approach which 
can be easily extended and implemented into 3D numerical models. 

We assume that a magnet is composed of perfectly aligned pointwise magnetic dipoles 
having the same magnetic moment. This assumption is well posed for modern manufac- 
tured permanent magnets, as follows from the final comparison between calculated and 
experimentally measured magnetic field. Take 2 as a direction of the unit magnetic dipolc 
m = (0,0, 1), then, a partial magnetic field in poin t r = (x, y, z) f rom a dipole located 
in point r' = (x',y',z') can be presented as, see e.g. Jackson ( 19991) : 



B'(r,r') = V x 

= V x 



m x 



R 

fl3 



= V x 



-m x V 





1 




mi 


(?) 


V x m 


= V 






R 







= V 



dz |r - 



(2.3) 



where R = r r' and R = r — r'|. Here we used R/R 3 — —V(l/R) along with few 
vector identities and omitted the constant /io/(47r). Then, the total field from a magnet 
occupying a space Q follows as: 

B(r)=/ B'(r,r')rfr' = v|-$(r), $(r) = / _*L (2.4) 
Jq dz J n |r - r'| 

The last integral can be computed analytically in some cases as we show below. For 
an arbitrary fi, the integration can be performed only numerically and then be once 
tabulated in a 3D array. This 3D array can be supplied into a numerical solver where a 
finite differentiation is applied to compute the external magnetic field. Another way is 
to differentiate analytically l/|r — r'| and then calculate numerically three integrals for 
each magnetic field component. 

When limits of the integration imposed by f2 are independent of each other, then 
the problem has an analytical solution by means of the indefinite integrals given in 
Appendix |A] 
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Figure 3. Comparison of computed (lines) and experimental (symbols) magnetic field intensities 
at different k = M y /L y , n = 0.02 (dot-dashed), 0.4 (dashed and symbols) and f.O (solid), 
2 = —0.7: B z component along x, y = 0, (a), B x component along x, y = 0, (b), B z component 
along y, x = 0, (c), B y component along y, x — 0, (d). 



As shown in Fig. f2J the magnetic dipoles are located in the region il = {\x'\ ^ 
M x , \y'\ ^ M y ,h ^ \z'\ ^ oo}, where 2h is a distance between north and south mag- 
netic poles. In the present paper, the condition \z'\ ^ h is assumed because the magnets 
used in the experiments are assembled onto a soft-iron yoke that closes magnetic field 
lines, i.e. the dipoles are effectively located from h up to infinity. By taking the corre- 
sponding derivatives of indefinite integrals given in AppendixjAj and after a few algebraic 
calculations one obtains: 



B x {x,y,z) 
B y (x,y,z) 
B z (x,y,z) 



1 

T 



arctanh 

k=±l j=±l»=±l 

— V V x ^ 

B 



V - j Mi 



E E 

k=±l j=±l»=±l 



(ijk) arctanh 

b~ E E E arctan 



r(i,j,k) 

x — iM-r 



fe=±l i=±l i=±l 



r(i, j, fc) 
'(af-*M s )(» 



A) 



(z - kh)r(i,j, k) 



(2.5) 
(2.6) 
(2.7) 



where r(i,j,k) — [(x — iM x ) 2 + (y — jM y ) 2 + {z — kh) 2 ] 1 ' 2 , and B Q is selected in such a 
way to have B z (0, 0, 0) = 1. Three-fold summation with the sign-alternating factor (ijk) 
reflects the fact that these equations are obtained by integrating according to Eq. (|2.4|) ■ 
Fig. [3] shows cuts of magnetic field intensities computed with equations (j2.5H2.7p for 
different k. Also, few experimental data (symbols) are presented for k — 0.4. There is 
a good agreement between experimental and analytical results. Moreover, one can see 
that the constrainment factor k affects mainly the spanwise distribution of magnetic 
field, whereas streamwise distribution changes only slightly. This is expected, because 
the length M x of the magnet is fixed, while the width M y = nL y is varied. It is necessary 
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to note that even for the broad magnet (n = 1.0) there is still a decline of B z near side 
walls, Fig. [3jc). 

Contour lines of the B z component in the middle plane (z — 0) at different k are also 
shown in Fig. [5Ja — c). 

Thus, we have demonstrated a self-consistent analytic approach to define arbitrary 
magnetic field configurations. It guarantees divergence- and curl- free requirements of 
B(r) and has a link with a clear physical model. 



3. Results and Discussion 

3.1. Stationary flow patterns in the middle plane 

In Section 13.11 we discuss characteristic stationary flow patterns which have been ex- 
tracted from 3D numerical results. Three-dimensionality of the simulations is of impor- 
tance to make these patterns stable, as we shall show later in Section l3~7l 

3.1.1. Streamlines for different constrainment n 

In our opinion the most striking effect of spanwise heterogeneity of the external mag- 
netic field is shown in Fig. 01 where flow streamlines in the middle plane, Fig. rjjd — /), 
are shown at the same flow parameters, N — 36 and Re — 196. To get an impression 
about the magnetic field configurations, the corresponding B z contour lines are shown in 
Fig. [4ja — c). Depending on the constrainment factor k, one observes the following sta- 
tionary flow patterns: a vortex dipole for the magnetic blade (k — 0.02), Fig. |4ja, d)> the 
stable six- vortex ensemble for the middle magnet (k = 0.40), Fig. \Mf>, e); and no vortex 
motion for the broad magnet (k = 1.0), Fig.|4jc, /). The projection of the magnetic pole 
onto the middle plane is shown by the bold solid line. 

Let us qualitatively explain these flow patterns. The case of the magnetic blade, 
Fig. [Ha, d), might be roughly understood by considering the limiting case of a Lorentz 
force which is pointwise in the spanwise direction. (By its definition, Fl = j x B, the 
Lorentz force is a volume force, however, if the distribution of heterogeneous magnetic 
field B is very sharp in space, the Lorentz for ce distribution is also very sharp.) As is 
known, see e.g. IVoropavev &: Afanasvev (1994), such an instant retarding force gener- 



ates vorticity which results in two counter-rotating vortices, a vortex dipole. Then, the 
induced vortices are advected and diffused downwards from the source of the force. In 
the case of MHD flow, these two vortices are fixed by a sidelong gradient of the mag- 
netic field, so they stay on the place . We should note th at the similar rec irculation has 
been ob tained numerically earlier bv lGelfgat et all ( 19781 ). and recently bv lCuevas et al. 



((200661) for a creeping flow. 

The case of the mid dle magnet, Fig. [4j&, e), is explained briefly in the recent work of 
Votvakov et al. ( 2007h by means of a mutual interaction of Lorentz and inertial forces. 



The first pair - inner magnetic vortices - is an inheritor of the vortex dipole as in 
Fig. E^d); the third pair - attached vortices - is of the same nature as a recirculation 
past a solid body, Fig. QJa), while the intermediate pair - connecting vortices - appear 
to make the coherent rotation of the magnetic and attached vortices possible, Fig.[lJ&), 
[4|e). It is clear why attached (hence, connecting) vortices are not induced in the case 
of the magnetic blade: it has a well streamlined shape, so there is no stagnancy region. 
This is in full analogy with the flow around a solid body where an appearance of the 
stagnancy region is strongly influenced by extent of streamlining. 

Finally, the case of the broad magnet, Fig. [4jc, /) shows no recirculation at given 
parameters, because the flow pattern is influenced now mainly by streamwise i nward and 
outward gradients of the magnetic field. As is known, see e.g. ISherclifj (|l962h . lKit etal 
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Figure 4. Contour lines of the transverse magnetic field component (a-c) and flow streamlines 
(d - /) in the middle plane (z = 0) at N = 36, Re = 196 and n = 0.02 (a,d), 0.4(6, e), and 
1.0(c, /). Magnetic pole is shown by bold lines. Dot-dashed(d), dashed(e) and solid(/) vertical 
cuts denote the location of velocity profiles of Fig. [5] 



the streamwise gradients and sidewalk of the duct are responsible 
all together for a Af-shaped profile of a streamwise velocity in the spanwise direction. 
The M-shaped profile can also develop a stagnant region in the middle of the duct 
at high interaction parameter N, however, no recirculation has been discovered until 
now. Moreover, as we shall prove in Section f3.4l no recirculation is possible if an external 
magnetic field is perfectly spanwise uniform. The present case of the broad magnet shows 
a slight decline of B z moving towards sidewalls, see Fig. [3Jc, solid line). In Fig. HJc), 
however, this decline is not enough to develop magnetic vortices at the given interaction 
parameter N = 36. We will see later that the critical interaction parameter N c , m which 
is needed to induce recirculation under the broad (k = 1.0) magnetic gap is more than 
one hundred. 
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Figure 5. Streamwise (u x , three lower curves, left y-axis) and spanwise (u y , three upper curves, 
right y-axis) velocities along spanwise cuts shown in Fig. \$^d ~ /). Middle plane, 2 = 0, and 
x = —4(a), x = 0(6), x — 4(c), x = 10(d). Constrainment factor k = 0.02 (dot-dashed), 
k — 0.40(dashed), and k = 1.0(solid lines), N = 36, Re = 196. Arrows show directions of 
spanwise redistribution of the flow. 

3.1.2. Velocity profiles at different n and the same N and Re. 

Fig. [5] shows quantitative data about stream- and spanwise velocity along different 
spanwise cuts. Different line styles correspond to the different k's. Arrows in Fig. [5] show 
directions of a spanwise redistribution of the flow. 

The largest braking Lorentz force, Fl = j x B, is generated at the front of a magnetic 
obstacle where induced electric currents j are maximum and the magnetic field B is 
strong. This results in a deformation of the incoming flow on the obstacle already at 
x = —4, Fig. [Sfa). The deformation consists of forming an inhomogeneous M -shaped 
profile of streamwise velocity u x and an appearance of a spanwise flow, i.e. u v , from the 
center to sidewalls, as shown by arrows in Fig. Ufa). 

Despite the order of magnitude difference in spanwise widths of magnets for k = 0.02 
(dot-dashed) and n — 0.4 (dashed), their ranges of reverse velocity (u x < 0, Fig. \Mp)) 
along the y-axis differ only by a factor one and half in the area of magnetic vortices. 
Positive streamwise velocities of the vortices are smoothly transformed into velocities of 
the external flow. Thus, for the wide range of k, the spanwise diameter of a single magnetic 
vortex remains nearly constant provided the vortex is not influenced by a sidewall. This is 
a manifestation of the fact that a decisive role for the vortex is the width of the spanwise 
decline of a magnetic field rather than the width of the magnet. 

By forming magnetic vortices, the M-shaped profile of streamwise velocity u x develops 
a negative value in the center, while the spanwise velocity u y changes its sign twice. There 
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is a redistribution of the flow from a vortex inner side to the centerline (y = 0), and from 
an external side of the vortex to the corresponding sidewall (y = L y ). On the contrary 
there are distinct zones past the magnetic vortices where a spanwise redistribution is 
absent (u y = 0), see \y \ ^0.8 for the magnetic blade (k = 0.02, dot-dashed) and \y\ ^ 1.5 
for the middle magnet (dashed), Fig. E^c). Also, the zone of u y = is observed behind 
attached vortices, \y\ < 1, Fig. \5^d), k — 0.4. 

Streamwise component of the velocity shows small changes with rising x, cf. u x in 
Fig. — d), except for a single peculiarity. For the case of a broad magnet (k = 1, solid 
lines), one observes that as x increases, the maximum of u x is shifted from sidewalls to 
the center due to diffusion of vorticity. In the same time, the u y component changes its 
sign in the region near sidewalls. 



3.2. Centerline profiles 

In this Section, taking as an example the middle magnet, k = 0.4, we show what forces 
are needed to induce vortices inside and past a magnetic obstacle. To reach this goal 
we shall analyze streamwise velocities along the centerline of the duct (y = z = Q), i.e. 
centerline profiles. 

The main difference between a magnetic obstacle and a solid body is the permeability 
of the obstacle depending on the retarding Lorentz force, Fl = j X B. This braking force 
is the largest in the center of the magnetic obstacle, where the magnetic field B reaches 
the highest intensity. The characteristic measure of the Lorentz force with respect to 
inertia is given by the interaction parameter N: higher the N, the stronger is Fl and 
less penetrable is the space under the magnets. Fig. [Sfa) is to illustrate this behavior: a 
dot-dashed line (N — 4) shows that the streamwise velocity is suppressed by approaching 
the magnet, however the braking force for the given N is not strong enough to reverse the 
flow. The reversion happens at higher N, see dashed (N = 9) and solid (N = 16) lines 
being negative inside the obstacle. The lowest velocity is marked in Fig.^a) by u Xi mi n , it 
becomes zero at a critical value iV CjTO , and at N N c ^ m one observes a recirculation - two 
inner magnetic vortices as shown in Fig. [4|g?, e). As the name points out, these vortices 
belong completely to the MHD flow, and have nothing common with a hydrodynamical 
flow around a cylinder. Rather these vortices are simi lar to those appearing under the 
action of a point braking force, see e.g. Afanasvev ( 20061 ). The concrete mechanism leading 
to the magnetic vortices is presented in Section 13.41 

Now, wc consider the behavior of centerline curves at fixed N and varying Re, Fig.^fe) 
and apply the analogy with ordinary hydrodynamics. Let us recall that the flow around 
a solid cylinder shows a stagnant region with two attached vortices when Re is slightly 
higher than a critical value. The typical centerline for this case is shown by inset of 
Fig. \§lb). The same situation can happen for the magnetic obstacle by increasing Re: 
the centerline curves past the magnetic gap become negative again as shown in Fig. [5^6). 
In fact, the attached vortices are induced past the magnetic obstacle, see Fig. [4je, third 
pair of vortices) analogous to those past a real solid body. 

Note, that the first minimum in a centerline profile is almost unperturbed when Re 
increases at fixed N, Fig. [5)^6), and this is another strong evidence that the vortices in- 
side and past the magnetic obstacle are of different physical origin. Since the interaction 
parameter is given as N = Ha 2 /Re, the magnetic vortices are enhanced by decreasing 
the flow rate (i.e. Re). On the contrary, as follows from Fig. [6^6), the attached vortices 
manifest themselves by increasing the flow rate provided that the intensity of the mag- 
netic field (i.e. Ha) as well as a spanwise magnetic field gradient have already enforced a 
reverse flow inside the magnetic obstacle. In other words, there is a qualitative distinction 
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(a), Re=100 (b), N=36 




Figure 6. Streamwise velocity varied along the centerline y = z = for the middle magnet 
(k = 0.4): effect of N at fixed Re = 100 (a); and effect of Re at fixed N = 36 (b). For (a): 
iV = 4(dot-dashed), 9(dashed), and 16 (solid lines); for (b): Re = lOO(dot-dashed), 144(dashed), 
and 196(solid lines). Vertical bold lines show borders of the magnetic gap, M x = 1.5. Inset in 
(b) is a centerline velocity profile for a flow around a circular cylinder with attached vortices. 
The streamwise velocity is normalized to the centerline velocity of unretarding flow. 



between magnetic and attached vortices: the former arise when Re decreases, and the 
latter arise when Re increases provided Ha is strong. 

It is easy to see in Fig. SJe) that magnetic and attached vortices are co-rotating in 
the same direction determined by the main flow movement. The only difference is the 
driving torque: this is the Lorentz force for the magnetic vortices, and the inertial force 
for the attached vortices. Because the magnetohydrodynamic and attached vortices are 
co-rotated, such a motion must be accompanied by a counter-rotation which produces 
the intermediate pair of connecting vortices, see Fig. HJe, second pair of vortices). The 
connecting vortices correspond to a local maximum on centerline curves, Fig. &b) behind 
the magnetic gap. 

Thus, N is responsible for the appearance of the magnetic vortices, while Re is respon- 
sible for the appearance of the attached vortices. The connecting vortices are necessary 
for the coherent rotation of the magnetic and attached vortices. 

Three decades ago Gelfgat et al. have attempted to reveal a kind of recirculatio n due 
to an external magnet ic field by both 2D numerica l simulation ( Gelfgat et al. ( 19781 )) and 
physical experiments ( Gelfgat fc Olshanskiil ( 19781 )). They saw a reverse flow numerically, 
and then designed a special experiment which did not confirm the recirculation. As follows 
from our results, the authors of the cited papers have not realized that they observed 
and discussed different phenomena in their numerical and experimental works. Their 2D 
numerical study neglected an inertial term, what corresponds to a creeping flow, and 
the observed reversion of the flow is just a sign of magnetic vortices. Recently, similar 
2D numerical work of ICuevas et all (|20066l ) showed the same eff ect in a creep i ng flow 
without side walls. Thus, the numerical work of Gelfgat precedes Cuevas et all ( 20066 ) 
by almost thirty years. 

The experimental report of Gelfgat &: Olshanskiil (|l978l ) contains a measured center- 
line profile (see Fig. 7b in the cited work), which falls behind the magnet, but does not 
approach negative velocities. In the cited experiment, parameters are N = 7.5 and 
Re = 3.73 x 10 5 , so the authors noticed that the observed falling velocity is due to 
inertial effects. Moreover, these experiments were performed for the turbulent flow be- 
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Figure 7. Plot (a): minimum of streamwise velocity on the centerlines curves (see Fig. |6ji) 
depending on the interaction parameter N = Ha 2 /Re obtained in simulations (closed symbols) 
and in experiments (open circles) for k = 0.4 and Ha — 140 (experiment), Re = 196 (squares) 
and Re = 100 (triangles). Plot (6): The critical interaction parameter N c , m depending on k for 
the appearance of magnetic vortices, 3d simulation, Re — 100. 



cause Re is rather high. To obtain the desired recirculation, the authors needed to just 
decrease further the inertial force by using a lower flow rate. This would have resulted in 
a higher N and, by exceeding some threshold, made impossible the turbulent pulsations. 
Then, the recirculation by the braking Lorentz force might have been revealed with the 
appearance of magnetic and attached vortices. Unfortunately, experiments with lower Re 
have not been performed, probably because the proper analysis of centerline curves (see 
Fig. 6) has been absent in those times. 

3.3. Existence regions of stationary flow patterns in parameter space 

A sign of the recirculation enforced by the Lorentz force is the negative value of tt x ,min, 
see Fig. [6f a) . We have performed a series of 3D simulations for the different interaction 
parameters N and constrainment factors k. An example of the series for the middle 
magnet (n = 0.4) is shown in Fig[7](a) . It is a ccompanied by e xperi mental data which 
were obtained by Oleg Andreev for the paper IVotvakov et all (|2007l ) but not reported 
there due to lack of space. Note, that the interaction parameter N — Ha 2 /Re is changed 
in experiments and numerics differently: experimentally by varying the flow rate (i.e. 
Re) at fixed external magnets (experimental Ha = 140), while numerically we kept Re 
constant and changed Ha. This is because of (i) experimental difficulties to work with a 
low flow rate in the closed channel, and (ii) numerical problems arising due to boundary 
layers at high Hartmann numbers. Moreover, a transition to turbulent regime is another 
factor that obscures stationary effects in numerics at high Re. 

As one can see in Fig. Eta), w x ,min is positive at low N, and then falls monotonically 
to a constant level showing perfect agreement between experiments and numerics for 
iV > 15. This excellent agreement confirms that the flow under magnets depends solely 
on N and magnetic field configuration rather than Re or inlet profile, provided, the 
external magnetic field is strong enough to sup press inward tu rbulent fluctuations. The 
latter condition was satisfied in experiments bv lAndreev et all (|2007l ). Also, a similarity 
of the electric potential map at equal N and different Re has been recently confirm ed by 
3D simulations and comparison with experiments, see Votvakov fc Zienickd ( 2007t ). The 
slight difference shown in Fig. Eta) at N < 15 is explained by the fact that the variation 
for N has been defined differently in experiments and numerics, as explained above. 

Dependence of w x ,min on N becomes negative at a critical value N c , m . This value 
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depends on the constrainment factor k. A series of 3D simulations has been carried out 
for the range 0.02 ^ n ^ 1 in the vicinity of u^ m - ul (N) ps and the results are given 
in Fig. which shows the separation line between regions for stable flow without 
and with magnetic vortice^j]. During this series of 3D simulations we did never find a 
hint that both flow patterns — the one with and the one without magnetic vortices — 
coexist for the same pair of parameters, i.e. both solutions are stable, but have different 
basins of attraction. However, we have checked different initial conditions only for some 
parameter combinations and did not carry out a systematic search with different initial 
conditions. Therefore, we can not conclude that the case of coexistence of two solutions 
is not possible. 

From the separation line of Fig. [7^6), as we have it numerically determined, the follow- 
ing trends are visible: lower the k, smaller the influence of sidewalls; the case of k — > 
corresponds to a free flow. Larger the k, more uniform is the braking Lorentz force in the 
spanwise direction. Therefore, in order to induce inner vortices at larger n it is necessary 
to apply a larger critical interaction parameter N c ^ m . For middle magnets, k^0.5, the 
critical value of N CjTn is of the same order of magnitude, N CtTn w 6, and for broad magnets, 
k ^ 0.8, it increases up to N c . m = 109 at k=1. The latter case has been earlier discussed 
in Fig. [4{c, /) for N = 36 as an example of a vortex-free flow pattern. We prove in next 
Section that any recirculation is impossible if the external magnetic field is perfectly 
spanwise uniform. 

3.4. The mechanism to induce recirculation. Vorticity, electric field, drop of pressure. 

One has to analyze the electric field inside the magnetic obstacle for different k in order to 
understand why the appearance of vortex motion is strongly dependent on the spanwise 
variation of the external magnetic field. 

Let us recall that the electric potential, (f>, is distributed according to the Poisson 
equation, A<fi = V • (u x B) = B ■ w, where w = V x u is the vorticity. Note that 
the projection of w on the externally fixed vector B plays a role of the induced electric 
charge density which creates the electrostatic field E = — V<j). It is possible to derive 
the direction of the field E from the knowledge of positive and negative extrema in the 
B • w distribution, by roughly assuming that these extrema can be approximated as point 
charges. As shown in electrostatics, the maximum (minimum) of rhs in Poisson equation, 
i.e. B z lu z , creates the minimum (maximum) in <f> distribution. 

For definiteness we consider the central spanwise cut, x — z = 0, where B x = B y = 
0, B z 0, and B • w = B z oj z « —B z du x /dy. This cut goes through the center of the 
magnetic gap and most expressively develops peculiarities typical of any spanwise cut. 
Vorticity uj z , the product of B z lu z , and resulting electric potential <f> along this central 
spanwise cut are shown in Fig. [8] for the broad (a) and middle (b) magnet. The magnets 
are depicted on the top and bottom by filled rectangles. 

The behavior of vorticity, uj z « —du x /dy, can be easily understood from the spanwise 
deformation of the streamwise velocity, i.e. M-shaped velocity profile, Fig. [5fc). Such 
a profile is characterized by two jets sidelong streamlining the magnetic obstacle, and 
each jet has the internal and external slope, hence, the vorticity alternates its sign by 

t It is tempting to call Fig. [7J&) a stability diagram. We renounce from doing this, because — 
strictly speaking — it is not proven that the change of topology of the flow patterns is caused 
by a change of stability. We find this assumption highly probable, but we did not compute 
eigenvalues of the Jacobi matrices for the stable flow patterns to check whether the transition 
is a bifurcation or not (see also our remark at the end of the introduction). We continue to call 
N c ,m the critical interaction parameter. This has to be understood not in the sense of a change 
of stability, but instead in the sense of separating stable solutions with a different topology. 
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Figure 8. Spanwise cuts, x = z = 0, of vorticity, u z (left axis, solid), magnetic field multiplied 
by vorticity, B z uj z (left axis, dashed), and electric potential, (j> (right axis), for k — l(plot a) and 

0. 4(plot b); N = 9, Re = 100. Filled rectangles at top and bottom show schematically magnets. 

passing the velocity maximum of the jet. The steepness of the external side of the jet is 
significantly higher than that of the internal side, especially for the larger k because the 
external slope is adjoined with no-slip sidewalls. Fig.[8]shows that the sidewalk vorticity, 

1. e. the external slope vorticity, is larger by an order of magnitude than the vorticity 
adjoining the center, i.e. the internal slope vorticity. To stress this fact, Fig. [8] shows 
both global and zoomed (insets) views for u> z and B z lo z . 

When the spanwise variation of the magnetic field B z (y) is weak, e.g. for the broad 
magnet (k = 1) shown in Fig. [5Ja), the product of B z uj z reflects completely the quanti- 
tative difference in lu z taken on the internal and external jet sides. Because the external 
w z is much higher than the internal u) z , the distribution of </> is determined mainly by 
the vorticity generated on sidewalls rather than by that produced inside the obstacle. 
Therefore, the negative (positive) external oj z induces strong positive (negative) poten- 
tial on the sidewalls, and <f>{y) drops monotonically along the y-axis, i.e. the spanwise 
electrostatic field E y = —dcj)/dy is always positive and does not change its sign between 
sidewalls. 

If the magnetic field is perfectly uniform along the y-axis, the electric potential distri- 
bution is completely governed by the sidewalls vorticity, and the alternating vorticity on 
the internal slopes of the M-shaped profile does not contribute to </>(y) even at large inter- 
action parameter N . As a result, the electrostatic field E ea (0, E y ,0) in the region of mag- 
netic obstacle is always directed opposite to electromotive force, u x B as (0, — u x B z , 0). 
By braking the flow, streamwise velocity u x and spanwise field E y accordingly approach 
zero by keeping their signs. A stagnant region can develop where streamwise velocity is 
small but still positive because the electric field is always in the same direction. Therefore, 
recirculation is impossible. 

Let us consider the case when the magnetic field intensity, B z (y) , drops from the center 
to sidewalls. Here, the high vorticity generated by sidewalls in the product of B z uj z is 
suppressed by a low intensity of the magnetic field. This is shown by the middle magnet 
(k = 0.4), Fig. Mj>)), the corresponding B z (y) is presented by dashed lines in Fig. [3jc). 
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The product of B z lo z resembles the u> z (y) behavior inside the magnetic obstacle, since 
there B z (y) ss const, i.e. the internal vorticity is well presented in B z lo z . Outside the 
obstacle, B z {y) rapidly decreases, so the external vorticity is remarkably weakened: as 
follows from Fig. [5^6), the dashed B z io z curve diverges from the solid io z curve. The largest 
difference between lj z and B z lo z are present on sidewalk where the external vorticity is 
very high due to no slip boundary conditions, while B z oj z at \y\ = L y is of the same 
magnitude because its alternating extrema at \y\ « M y correspond with borders of the 
magnet. Therefore, the dependence of 4>(y) is allowed to be influenced by the internal 
slopes of the M-shaped profile, which produces extrema on <fi(y) at \y\ « M y alternating 
with the sign of </> on sidewalls. These internal extrema on </>(y) reverse the spanwise 
electric field i.e. E y = —d<fi/dy becomes negative inside a finite magnetic obstacle. 

The reverse spanwise electrical field appearing for spanwise decaying magnetic field is a 
necessary but not sufficient condition for the appearance of reverse flow. It is a necessary 
condition, because otherwise no Lorentz force pointing in negative x-direction would 
appear. This becomes clear looking at equations (|2.ip and (|2.2|) for Lorentz force and 
current density: j xB — — V</>xB — £> 2 u+(u-B)B. Considering again the central spanwise 
cut, x — z — 0, where B = (0,0, B z ), u w (u x ,0, 0) and j « (0,j y ,0), one can write the 
x-component of the Lorentz force as follows: (j x B)^ rj j y B z = E y B z — B 2 z u x . The third 
term has vanished, the second term corresponds to a factional force proportional to the 
actual velocity. The only term, which can provide a driving force for the reverse flow is the 
first term, and this only then when E y is negative, i.e. when a reverse spanwise electrical 
field exists. Only the existence of a reverse spanwise electrical field is not sufficient to drive 
reverse flow, because it has to reach a minimal strength. This becomes clear from the 
y-component of Ohm's law for the electrical current density at the central spanwise cut: 
j y = E y — u x B z . To see this, one needs an additional information concerning the global 
behavior of the electrical current in our system. The electrical current is organized in two 
horizontal loops, with a strong negative spanwise current (i.e. j y < 0) directly under the 
magnet, which divides at the sidewalls into one loop before the magnet and the other 
behind the magnet (see figures 16, 17 and 18 of our former paper IVotvakov fc Zienickd 
(2007)). The sign of j y does not change, when the interaction parameter is increased. 
When E y just changes sign at a certain value of N at the start of electric field reversal, 
—u x B z has to be negative. This is only possible for (still) positive u Xl and consequently 
N has to be further increased until the reverse spanwise electrical field becomes so strong 
that it equals the value of j y . This characterizes, in fact, the critical iV Cjm , for which the 
velocity u x is exactly zero. Let us denote this value of the reverse spanwise electrical 
field by E VtC . A further increase of the interaction parameter makes E y stronger negative 
than j v with the consequence that u x has to be smaller than zero, which is equivalent 
to the existence of reverse flow and recirculation. Thus, summing all up, the role of 
spanwise magnetic field gradient is to suppress the external vorticity and promote the 
internal vorticity in the product of B ■ u>. When by a further increase of the interaction 
parameter the internal spanwise electrostatic field is reversed stronger than the critical 
Ey tC , the recirculation appears. 

Fig. (9](a) shows E y = —d(f>/dy for the different magnets taken in the center of the 
magnetic gap, x — y = z = 0, as a function of the interaction parameter N normalized 
with the critical value N c>m . Letting the appearance of recirculation be the referred flow 
pattern, the ratio N/N CyTn characterizes the power of recirculation for various k. One 
can see that the curves are close to each other when N/N c , m > 1. Moreover, by taking 
N/N Cym = 1 one finds that the critical E y<c w —0.17 is independent of n. That is, the 
critical magnitude of the spanwise electrostatic field, E y _ c is a universal parameter which 
is the same for different magnetic field configurations. This fact is clear because B z = 1 
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Figure 9. For various magnetic field configurations, k = 0.4 (open circles), 0.6 (filled circles), 
0.8 (open squares), 1.0 (filled squares), the dependencies on N/N c , m in the center of magnetic 
obstacle, x = y = z = Qoi the spanwise electrostatic field, E y = ~dcj>/dy, (plot a), and the 
normalized resistance to the flow, r = (~dp/dx) / N c , m , (plot b). 



and j y — Ey_ c when the recirculation starts, i.e. u x = u y = 0. So, E VtC is indeed the critical 
braking Lorentz force because Fl x — j y B z — E VtC . 

Let r — — (dp / dx) / N c , m atx = y = z = 0be the resistance to the flow inside the magnetic 
obstacle. At the beginning of recirculation the Navier-Stokes equation is simplified up 
to dp/dx + N c ^ m E V:C — 0, so r c = —E VtC . The behavior of r as a function of N/N C:m 
for various n is shown in Fig[!J&). One can see that the change of the flow regime is 
accompanied by the change of the slope in the resistance oir(N/N Cjm ). For those k where 
the sidewall influence is insignificant, the appearance of recirculation at N/N CjTn > 1 
results in a drop of the resistance despite the fact that N increases. 

The second effect of a drop in the resistance due to magnetic recirculation has the same 
explanation as the drag crisis well known in ordinary hydrodynamics. It appears in the 
wake of a circular cylinder when the boundary layer on the cylinder surface undergoes a 
transition from the laminar to turbulent mode. It causes a substantial reduction in the 
drag force analogous to that in the wake of the magnetic obstacle. 



3.5. Similarity of MHD flows inside the magnetic obstacle 

As follows from Eq. (|2.ip the viscous force Au is scaled by Reynolds number Re, therefore 
at high Re and far from the walls it plays a minor role. As a result the flow must be 
governed by the interaction parameter N defin ing the ratio between magne tic and inertial 
forces. This peculiarity was marked already in lVotvakov Sz Zienickd ([20071 ) by comparing 
experimental and numerical electric potential distributions found at different Re and 
similar N. 

In the present paper, Fig. 1101 demonstrates the same behavior for the broad (curves 1 
and 3) and middle (curves 2) magnet. Solid and dashed curves 1, 2 have the same N but 
Re = 100 for solid and Re — 400 for dashed lines. 

It is expected after discussing Fig. [5] that the results given in Fig. [TU] are similar 
even for different constrainment factors k and interaction parameters N provided the 
corresponding ratios N/N c ^ m are equal. Curve 2 for the middle magnet is plotted for 
N = 9, i.e. N/N Cym (n — 0.4) = 9/6.8 = 1.32 and curve 3 for the broad magnet - for 
N = 144, i.e. N/N c , m (K = 1.0) = 144/109 = 1.32. One can see that in the central part 
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Figure 10. Spanwise dependence [x = z = 0) of electric potential (plot a) and streamwise 
velocity (plot b) for the broad (n = 1.0, N Cm = 109, solid and dashed curves 1,3) and middle 
(k — 0.4, N Cm = 6.8, solid and dashed curves 2) magnetic obstacle. N = 9 (curves 1,2) and 
144 (curves 3), Re = 100 (solid) and 400 (dashed lines). Notice that for the curves 2 and 3 the 
constrainment factors k are different, while the relations N/N C:m « 1.32 are equal. 



\y\ < 2 being approximately the spanwise dimension of the less magnet (M y = 0.4x 5 = 2) 
both electric potential (a) and streamwise velocity (b) are also close to each other, cf. 
curves 2 and 3. 

3.6. 3D peculiarities of the flow 

Fig. [TT1 shows different 3D space perspective views of the six- fold vortex pattern for the 
middle magnet (k = 0.4), N = 16, and Re = 196. For easy visualization, plot a shows the 
XY projection of all vortices, i.e. magnetic, connecting, and attached, while other plots 
present only part of the recirculation. Plot b shows the XZ projection for magnetic and 
attached vortices, plots (c — d) are both YZ projections of the magnetic, and attached 
vortices, correspondingly. 

The following 3D features are observed in Fig. [IT] (i) a helical motion inside every 
vortex; (ii) the axes of rotation of the magnetic vortices are parallel to the lines of 
external magnetic field; (iii) the two magnetic vortices are closely adjoined to each other 
and taken together they form a barrel extending in its central part and located mainly 
between magnetic poles; (iv) two helices of attached vortices are not adjoined and are 
arched along x-direction. 

Additionally, Fig. [12] present flow streamlines constructed from u y and u z components 
of velocity field for few characteristic vertical slices: at front of the magnet (x = —4, 
plot a), and the vertical cuts of magnetic (ie = 0, plot 6), connecting (a; = 4, plot c) and 
attached (x = 8, plot d) vortices. These (u v , u z ) streamlines are not real moving lines 
of fluid particles in the vertical slices, because the particles have also the u x velocity 
component. Instead, these streamlines are to show ascending and descending paths of 
fluid particles in vertical slices as projection. 

At x = —4, Fig. [TSb . one observes the braking effect of the Lorentz force making the 
flow streamlines diverge from the central point y = z = 0, where the greatest change 
of du x /dx = —(duy/dy + du z /dz) occurs. Because the u x component is not involved in 
the plot, this effect looks as the source of the flow. Also there is a motion in z direction 
caused by a tendency to form Hartmann layers, see Votvakov fc Zienickel ( 2007 ). 



At x = 0, Fig. [T2b . one firstly notices the abrupt change of the motion in y direction 
at | y | w 2. This helps us to see the vertical borders of the magnetic vortices where the 
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(b), XZ view 




(c), YZ view 




Figure 11. Perspective projection of 3D vortex disposition for middle magnet and N = 16, 
Re = 196: XY view, magnetic, connecting and attached vortices (a), XZ view, magnetic and 
attached vortices (b), YZ view, magnetic vortices (c), YZ view, attached vortices (d). 



(a), x=-4, N=16, Re=196 



(c), x=4, N=16, Re=196 





(b), x=0, N=16, Re=196 



(d), x=8, N=16, Re=196 





Figure 12. Flow streamlines constructed from u y , u z velocity components in vertical slices 
x = -4(a), 0(6), 4(c), and 8(d); iV = 16, Re = 196. 



20 E. V. Votyakov, E. Zienicke and Yu. Kolesnikov 

u y reverse its sign. Moreover, a barrel shape of the magnetic vortex dipole can be clearly 
seen. Directions of arrows on the flow streamline demonstrate a drift of the flow from 
top/bottom walls towards the center inside the magnetic vortices, and a drift in the 
opposite direction outside. Such a drift is a reflection of helical motion which can be seen 
in perspective in Fig. [Til 

Similar to Fig. 112b . sharp borders between the main flow and connecting (x = 4, 
Fig. fT2b ) and attached (x — 8, Fig. \V2ti ) vortices are observed. One also observes a 
remarkable vertical drift to and from the top and bottom walls depending on the direction 
of horizontal rotation. The vertical drift at x = 4 and x — 8 is more pronounced than that 
of at x = because the magnetic field has a much larger intensity in the case of x = 0. 
The new phenomena compared to Fig. 112b are swirls in the corners of the duct shown in 
plots (c, d). Similar swirls can also appear when a MHD flow has no recirculation. They 
arise due to the inertial f orce and destruction of Hartm ann layers resulting from a decline 
of the magnetic field, see Votvakov k. Zienicke] ( 2007 ). 



The vertical drifts shown in Fig. I12t > — d are caused by rotation of the fluid and can 
be unders t ood a s a manifestation of the hydrodynamic Ekman pumping effect, see e.g. 



Davidson! (|200lh . There are six rotating columns in the MHD flow studied here, and 
each one has its own primary horizontal motion and secondary vertical drift. Taken all 
together it gives rise to the helical motion clearly observed in Fig. [IT] 

The whole 3D space trajectory for the infinitesimally small volume of the fluid - a fluid 
particle - can now be described in the following way. Far upstream from the magnetic 
system, the particle moves straight under the pressure gradient. Approaching the region 
of influence of the magnetic obstacle, the particle turns aside towards the closest corner 
due to the action of the braking Lorentz force and reaches a boundary layer. If the particle 
is not captured by recirculation, it then passes the region of the magnetic obstacle in 
the bulk of two jets sidealong streamlined the obstacle. If the particle is captured by 
recirculation, it first goes down (up) from the top (bottom) towards the middle plane 
in the helix of the magnetic vortex. In the middle plane, this trajectory is taken by the 
helix of the closest connecting vortex and passes helically towards the top (bottom) wall 
to come again to the boundary layer and dissipate the kinetic energy. Then, the particle 
can be caught by the helix of the attached vortex and be drifted slowly to the middle 
plane where it finally becomes free to go downstream from the magnet. 

Recirculation under magnetic poles brings new details into top and bottom boundary 
layers perpendicular to the magnetic field. These layers are Hartmann layers in the case 
of a constant magnetic field, Ekman layers in the case of an axisymmetric rotating flow 
bounded by a fixed horizontal plate, and Ekman-Hartmann lay ers when the axisymetri c 
rotating flow is s ubjec t to a constant vertical magnetic field, lAcheson fc Hide! (|l973h . 



Desiardins et al. ( 1999^3- We checked velocities profiles u x {z) in the region of magnetic 



vortices and found no satisfactory agreement with the current theory for Ekman-Hartman 
rotating flows. The reasons for this disagreement are probably that the recirculation 
induced by a heterogenous external magnetic field shows details which are not compatible 
with the present analytic theory: (i) the rotating flow is not axisymmetric, (ii) neither 
magnetic field nor angular rotation are constant, (iii) this is a system of six rotating flows. 
The quantitative analysis of these layers requires further detailed investigation. At the 
moment we can conclude only that these layers are important to stabilize recirculation 
even at high Reynolds numbers as shown in next Section. 



f We thank the Referee who brought these works to our attention. 
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3.7. Recirculation in a 3d flow versus a vortex generation in a 2D flow 

This paper is devoted to a stationary 3D MHD flow. Recently, ICuevas et all ( 2006a ) 
have reported a vortex generation in a 2D flow induced by a magnetic obstacle. In this 
Section, we discuss how their 2D results are related to our 3D results. 

In our opinion, it is an open question on how to build a 2D model at high Re. Two- 
dimensionality assumes that the flow rate is kept constant in the plane under consid- 
eration. This assumption is certainly wrong in the case of a local magnetic field, where 
Hartmann layers are formed (destroyed) under inward (outward) magnetic field gradient, 
i.e. streamwise velocity profile in the transverse direction is becoming more (less) flat, 
th erefore the fluid must go out of (go into) the plane, see discussion about Fig. 8, 9 given 
bv lVotvakov fc Zienickd <|2007h . 

On one hand, the friction imposed by a no-slip wall stabilizes a flow because it provides 
a sink for kinetic energy. For instance, in ordinary hydrodynamics there are numerous 
experiments which ill ustrate that t he co n finement of the endplate s increases the stability 
of the wake, see e.g. IShair et all (|l963l ). iNishioka &: Sato! (1974]), Gerich fc Eckelmannl 
(|!982l ). lLee fc BudwiJ (|199lh . These examples include the delay of the critical Reynolds 
number for vortex shedding and the extension of the Reynolds number range for a 2D 
laminar shedding. 

On the other hand, the case of a MHD flow under a local magnetic field always requires 
to take into account top and bottom endplates because these plates carry magnetic poles. 
In reality, one can increase the distance between poles to have more two-dimensionality 
in the middle plane, but this automatically decreases the degree of space heterogeneity of 
the external magnetic field. So, it becomes an issue whether it is practically possible to 
design a strongly heterogenous magnetic field having a large distance between magnetic 
poles. 

There is a method to take no-slip top /bottom plates into consideration by means of 



the Hartman n friction term , see e.g . Lavrentiev et al. ( 1990Tj . This assumption has been 



used also bv ICuevas et al. ( 2006a ). However, this averaged approach is well validated 
only in the case of a small flow rate, i.e. low Re, even when the magnetic field is not 
strongly varying. It follows from the fact that past the magnetic obstacle the Hartmann 
friction takes a form of the Hclc-Shaw friction based on the assumption that velocity 
is parabolic along the z-axis. This is well validate d only f or vis cous flows where the 
vorticity advection does not play a significant role, iRieg e3 (|l938l ). The formal vertical 
Hele-Shaw friction term inserted into the 2D Navier-Stokes equations does not describe 
quantitatively the behavior of the system with high Re. We believe that this friction has 
no significant influence when the advection is strong. 

Nevertheless, one can consider mathematically what happens in a 2D free (no side 
walls) MHD flow where the flow rate is kept constant. This mathematical problem has 
been addressed bv lCuevas et all ( 2006al ). 

The cited 2D MHD flow had no walls, and, hence, no sinks for kinetic energy except an 
internal viscosity which is negligible when Re is high. In such a flow, any vortex pattern 
cannot be stabilized, and one would not observe any stationary recirculation. Instead 
one obtains the nonsteady vortex generation, hence, it is worth to discuss a developing 
flow, i.e. the flow starting from the constant velocity field u(r)| t=0 = (u x , 0) . We ha ve 
perform ed a few runs of 2D simulation in order to reproduce results of ICuevas et 



( 2006al ) and find out whether the six- fold vortex pattern shown in Fig. QJ6) is of general 
matter. These results are presented in Fig. Q2] 

One observes, at the initial times, subsequently symmetric magnetic vortices (Fig. ll3f a)). 
and then a six-vortex pattern (Fig. [T37 &. c)). This pattern quickly grows in size due to 
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(a), t=1 (d), t=90 




Figure 13. Flow streamlines in a 2D nonsteady free flow: t = 1(a), 10(6), 20(c), 90(d), 120(e), 
150(/); u(r)| t=0 = (1.5,0), N = 30, i?e = 100. 
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Figure 14. Vorticity field corresponding to Fig. I13f f). Dashed lines are plotted for negative, 
and solid lines are plotted for positive vorticity. 



inertia, (Fig. [T3]) (&. c, d), in such a way that attached vortices gradually swell and reach 
a dimension much larger than that of magnetic and connecting vortices. When attached 
vortices exceed their critical size, they become unstable and lose symmetry, Fig. [l~37 e). 
This gives rise to the Karman vortex street, Fig.[13]J/), illustrated also by a vorticity con- 
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tour plo t in Fig.Qj] The Karman vortex street has been reported already by lCuevas et 



(|2006a[ ). while the preceding temporal evolution, Fig. [TST a — /), is explained here for the 
first time. 

ft is known that the attached vortices in the flow past a solid body can be discovered at 
any Reynolds number Re at the initial instance of time before vortex shedding starts. One 
can see that the same situation takes place in an initial How past the magnetic obstacle. 
The difference is that instead of just attached vortices one can observe a six-fold vortex 
pattern developing from the recirculation under the magnetic gap. 



4. Conclusions 

We have reported the results of a 3D numerical study about a stationary liquid metal 
flow in a rectangular duct under the influence of an external magnetic field. The inter- 
action parameter N, Reynolds number Re as well as magnetic field configuration have 
been systematically varied. Whenever it has been possible, the numerical results have 
been quantitatively compared with experimental ones. 

First, an analytical, physically consistent and simple model for the external magnetic 
field has been derived. Parameters of the model are geometric dimensions of a region 
occupied by external magnets. This model has been successfully verified by comparing 
with experimentally measured data and then used through the paper by varying the 
constrainment factor, n = M y /L y , the ratio between spanwise dimension of the magnet, 
2M V , and width of the duct, 2L y . 

One can classify the following three typical flow structures into which a stationary 
MHD flow is organized, depending on the interaction parameter N as well as spanwise 
magnetic field heterogeneity. 

The first structure, attributed to low degree of field heterogeneity, is characterized by 
the significant electromotive force which opposes the electrostatic field. The characteristic 
pattern for this case is the Hartmann flow. If magnetic field is uniform, this regime takes 
place at all N, otherwise it is realized at such N when the field heterogeneity is not 
strong manifested. 

The second stationary structure is perfectly observed for the fringing magnetic field, 
i.e. when the intensity of transverse magnetic field is varied slowly in the spanwise and 
strongly in the streamwise direction. The flow pattern is given by the M-shaped stream- 
wise velocity profile without recirculation inside the magnetic gap. Here, the electromo- 
tive force and the electrostatic field can either be opposed or be in the same direction, but 
the direction of the electromotive force is always in the direction of the electric current. 
The loops of the electric current are located mainly in the horizontal plane. 

The decisive condition for the appearance of the third flow structure is a strong span- 
wise variation of the magnetic field which induces recirculation inside the magnetic obsta- 
cle. The recirculation starts when the reverse electrostatic field prevails a critical value. 
Here, the electromotive force is opposed to the electrostatic field and the direction of 
the electric current. For this recirculation regime, the intensities of the reverse flow, ob- 
tained in 3D numerically and by physical experiments have been compared, and a good 
agreement has been observed. 

The existence regions of stable stationary flow patterns, that is the dependence of 
the critical interaction parameter, N Cjm , to induce recirculation on the constrainment 
factor, k, has been calculated and discussed. It has been made clear that no recirculation 
is possible for perfectly spanwise uniform external magnetic field. Moreover, the MHD 
flows for various k have been shown to be similar provided they are of the same ratio 
N/N c , m . 
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Finally, 3D features of the flow under consideration have been discussed and it has been 
demonstrated that the magnetic vortices are stable in their disposition. This is contrary 
to a 2D numerical study where the stationary recirculation is possible only in a creeping 
flow while at higher flow rate the recirculation develops a vortex shedding. Nevertheless, 
one can see all these vortex patterns in a 2D nonsteady flow at initial times. 

We do not want to close this work without some hypotheses on the nature of the tran- 
sitions that we have found in this system, namely, (1) the transition between streamlining 
flow and the flow with magnetic vortices when N c , m and k are varied, and (2) the tran- 
sition from the two-vortex pattern (only magnetic vortices) to the six-vortex pattern 
(magnetic, connecting and attached vortices) when the Reynolds number is varied. The 
first transition has a high probability to be a topological change of the same solution, 
which is stable in the whole space of initial conditions. Increase of the strength as well 
as the increase or decrease of the spanwise inhomogeneity of the magnetic field seem 
us to be topological changes of the force field resulting in a topological change of the 
stationary stable solution. This would be consistent with the fact, that we never found 
the different flow patterns coexisting for the same parameter pair (N c . mi Hi). The second 
transition mentioned is to our opinion analogous to the appearance of attached vortices 
behind a solid obstacle as known from usual hydrodynamics. The magnetic vortices act 
as an obstacle for the flow. For increasing Reynolds number a shear instability arises 
resulting into the formation of attached vortices, which for consistence of the flow have 
to be accompanied by the connecting vortices. These impressions on the questions of 
stability that we got from our research on the considered system, nevertheless, remain 
to be proved rigorously. 

The authors express their gratitude to the Deutsche Forschungsgemeinschaft for finan- 
cial support in the frame of the " Research Group Magnetofluiddynamics" at the Ilmenau 
University of Technology under grant ZI 667. The simulations were carried out on a JUMP 
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at the Forschungszentrum Jiilich. We are grateful for many fruitful discussions with An- 
dre Thess. Special thanks go to our experimental colleague Oleg Andreev, for an always 
close exchange of thoughts and for providing the experimental data to compare with our 
numerical results. 



Appendix A. Indefinite integrals 

We take the following notations: r = (x,y,z), r' = (x',y',z'), Ar = \r — r'\ — 

1 /2 

[Ax 2 + Ay 2 + Az 2 ] 1 and Ax = x - x' , Ay = y - y', Az = z - z'. Then, the indefinite 
integration over z' gives 




(Al) 



the indefinite integration over z' , y 1 gives 




+ Ayarctanh 



Az 



+ Az arctanh 



Ay 



(A 2) 



Ar 



Ar 



and finally the indefinite integration over z',y',x' gives 
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&*>y^(x,y,z,x',y',z') 



1 



dx'dy'dz' 



Ax 2 arctan 



H — Ay arctan 
2 y 

1 A 2 

+ 2 arctan 



AyAz 
AxAr 
AxAz 



AyAr 
AxAy 



AzAr 



AyAz arctanh 
AxAz arctanh 
AxAy arctanh 



Ax 



Ar 
Ay 



Ar 
Az 
Ar 



(A3) 



To confirm, one can use a commercial program, e.g. Mathematica, to analytically 
differentiate (|A 3[) back and get (|A 1[) finally. The original way to calculate (|A 1[) - (|A 3p is 
quite cumbersome, and is not provided here. 
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